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ABSTRACT 
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^ . Using early data from the Infrared Array Camera (IRAC) on the Spitzer Space 

^ ! Telescope, taken for the Great Observatories Origins Deep Survey (GOODS), we 
identify and study objects that are well-detected at 3.6/im, but are very faint 

. (and in some cases, invisible) in the Hubble Ultra Deep Field (HUDF) ACS and 

O ! NICMOS images and in very deep VLT i^s-band imaging. We select a sample of 

! 17 objects with fui^>-^l^i^) / fviz^^o) > 20. The analysis of their spectral energy- 

OO 
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Tj- . cannot be satisfactorily explained without a well-evolved stellar population. We 
O 

:^ 
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find that most of them can be well fitted by a simple two-component model, 
where the primary component represents a massive, old population that domi- 
nates the strong IR emission, while the secondary component represents a low 
amplitude, on-going star-formation process that accounts for the weak optical 
^ ! fluxes. Their estimated photometric redshifts (zp) range from 1.6 to 2.9, with the 

median at Zp = 2.4. For the simple star formation histories considered here, their 
^ I corresponding stellar masses range from 0.1 to 1.6 x IO^^Mq for a Chabrier Initial 

^ I Mass Function (IMF). Their median rest-frame i^^-band absolute magnitude is 

—22.9 mag in AB system, or 1.5 x L*{K) for present-day elliptical galaxies. In 
the scenario of pure luminosity evolution, such objects may be direct progeni- 
tors for at least 14 to 51% of the local population of early-type galaxies. Due 
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to the small cosmic volume of the HUDF, however, this simple estimate could 
be affected by other effects such as cosmic variance and the strong clustering of 
massive galaxies. A full analysis of the entire GOODS area is now underway to 
assess such effects. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: 
luminosity function, mass function — infrared: galaxies 

1. Introduction 

The Spitzer Space Telescope^ the fourth and last of NASA's Great Observatories, pro- 
vides order-of-magnitude improvements in capabilities of studying the infrared sky over a 
wide wavelength coverage (Werner et al. 2004). IRAG, one of the three instruments on- 
board Spitzer, is a four-channel camera that simultaneously takes images at 3.6, 4.5, 5.8, 
and 8.0/im (Fazio et al. 2004). Extremely deep observations with IRAG are a key component 
of the GOODS Spitzer Legacy program. By the time of writing, the first epoch of GOODS 
IRAG observations has been finished. With a typical exposure time of 23 hours per pixel, 
these data have already imaged the 3.6-8.0/im sky to an unprecedented depth. About 1/3 
of the GOODS area, including the Hubble Ultra Deep Field (HUDF ^; PI. S. Beckwith), has 
been observed in all four IRAG channels. 

A subject of immediate interest is whether these IRAG data have revealed any objects 
with unusual properties at these largely unexplored depths and wavelengths. In this paper, 
we discuss a particular population that are bright in all IRAC channels, but arc extremely 
faint or even invisible at optical wavelengths. The red colors of these objects are reminiscent 
of those for the so-called "Extremely Red Objects" (EROs; e.g., Elston ct al. 1988; McCarthy 
et al. 1992; Hu & Ridgway 1994; Thompson et al. 1999; Yan et al. 2000; Scodeggio & 
Silva 2000; Daddi et al. 2000), which are commonly selected based on R — K or I — K 
photometry. We will therefore refer to our objects as IRAG-selected Extremely Red Objects 
(lEROs), and will discuss their possible connection to conventional EROs in this paper. To 
better constrain their optical fluxes, we concentrate our discussion on the area defined by 
the HUDF, where the deepest optical data are available. Throughout this paper, we adopt 
the following cosmological parameters: VLm = 0.27, Vl^ = 0.73, and Hq = 71 kms~^ Mpc~^. 
All magnitudes are in the AB system unless specified otherwise. 



^See http://www.stsci.edu/hst/udf 
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2. Data, Photometry and Sample Definition 

The data from the first epoch of GOODS IRAC observations on the CDF-S are described 
by Dickinson et al. (2004, in preparation). The nominal exposure time per pixel is about 
23.18 hours in each channel. The southern 2/3 of the entire field has been covered by 3.6/im 
and 5.8//m channels, while the northern 2/3 has been covered by 4.5/xm and S.O/im channels. 
The middle 1/3 of the field, which includes the HUDF, has been observed in all four IRAC 
channels. The images were "drizzle" -combined (Pruchter & Hook 2002), and the pixel scale 
of the final mosaics is 0.6", or approximately half of the native IRAC pixel size. For an 
isolated point source, the formal detection limits {S/N = 5) from background shot noise 
in regions with full exposure time range from 0.11 (in 3.6yum) to 1.66 /iJy (in S.O/im). In 
practice, crowding and confusion influence the detection limits. As discussed later, we will 
restrict ourselves to objects which are reasonably well isolated. 

We used SExtr actor (Berlin & Arnouts 1996) in double-image mode to perform matched- 
aperture photometry. We detected objects in a weighted average of the 3.6/im and 4.5/im 
images in order to provide a single catalog that covers the entire GOODS area with reason- 
able uniformity. A 5 x 5 pixel Gaussian convolving kernel with a FWHM of 1.8" was used 
for detection. We required a detected object have a minimum connected area of two pixels 
(in the convolved image) that were 1.5 cr above the background. We adopted the photo- 
metric calibration constants provided in the image headers generated by the Spitzer Science 
Center IRAC data processing pipeline. The "MAG_AUTO" option was used throughout. 
Photometric errors were estimated using realistic noise maps generated as part of the data 
reduction process, but these include background photon noise only. Crowding and source 
blending will generally increase the photometric uncertainties, although in this paper we 
will limit ourselves to reasonably well isolated objects. In total, there are 552 IRAC sources 
detected within the solid angle covered by the ACS images of the HUDF (10.34 arcmin^ in 
size) . 

We used the HUDF ^gsQ-band-based catalog of Yan & Windhorst (2004) for optical 
identification. The magnitudes in this catalog are matched- aperture "MAG_AUTO" mag- 
nitudes. The source matching was done by identifying the closest ACS object within a 1" 
radius from the IRAC source centroid. The S.dfjLm-to-Zsso flux density ratio (r(3.6/2;) = 
fui^-^jJ-m) / f,^{zs5o)) histogram of the matched sources has its peak at 1.7. We define lEROs 
as objects that have r(3.6/z) > 20. The unmatched IRAC sources automatically satisfy this 
criterion. No attempt has been made to correct the differences in flux measurement that 
caused by the different apertures used in generating the IRAC catalog and the ^gsg-band 
catalog, since such differences are at most ~ 0.1 mag (based on Monte Carlo simulations 
using artificial sources and tests using larger photometry apertures) and thus will not affect 
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the selection for these extreme objects. 

The major difficulty for reliable lERO selection and photometry is contamination due to 
source blending. Here, we adopt a conservative approach, by visually inspecting each source 
in both the IRAC and the ACS + NICMOS HUDF images, and considering only those 
sufficiently isolated objects that crowding should not significantly affect the photometry. We 
inspected 75 IRAC sources that met the criterion T{'i.Q/z) > 20, and rejected 58 that had 
nearby objects that might have influenced the IRAC photometry, or that were clearly the 
blended emission from two or more sources as revealed by the ACS images. The remaining 
17 objects constitute our final sample. Fifteen of these objects do not have any companion 
within a 1" radius from the IRAC centroids as seen in the ACS images. The other two have 
one or two close neighbors in the ACS images. However, we believe that their identifications 
are secure for two reasons: in the ACS bands these neighbors are all at least 1.0 mag 
fainter than the identified counterparts, and are all centered more than 0.5" from the lERO 
centroids. 

Table 1 gives the coordinates and photometric properties of these 17 lEROs. The source 
fiuxes are in the range 0.64 < /,^(3.6//m)///Jy < 11. Among these objects, 12 are within 
the field coverage of the HUDF NICMOS Treasury Program (Thompson et al. 2004, in 

preparation), of which 11 (including the two objects that have no optical counterparts) arc 
identified in both the Jno and the i/ieo bands. The unidentified one is the faintest IRAC 
source in our sample (object #10). In addition, we also searched for their counterparts in the 
deep f^s-band images obtained by ISAAC at VLT, which were taken as part of the ground- 
based supporting data for GOODS (Giavalisco et al. 2004). Fifteen of the lEROs have 
counterparts in the Kg images, and the other two are not seen because they are too faint. 
The Jiio, HiQo and Kg magnitudes of the lEROs are also listed in Table 1 when available. 
As an example, the image cut-outs of one of these objects are displayed in Figure 1. 



3. The Nature of lEROs 

Without exception, the lEROs that have optical counterparts all show a monotonic 
trend of increasing in fiux from i775-band to 3.6//m-band, and most of them show this trend 
starting in the S435-band. This fiux-increasing trend closely resembles that of the "unusual 
infrared object" found by Dickinson et al. (2000) in the HDF-North, whose nature remains 
uncertain. 

Here, we look for the simplest explanation for the lERO population as a whole. The 
red nature of these objects immediately leads to a number of broad possibilities: Galactic 
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brown dwarfs, old and/or dusty galaxies, and objects at extremely high redshifts. We can 
rule out the possibility that they are stars in the Galaxy, because they are all resolved in the 
HST images (ACS and/or NICMOS). Furthermore, as 15 out of the 17 objects are visible in 
the ultra-deep optical images to at least ijj^-hand, we can also reject the very-high-z {z > 7) 
interpretation. In fact, most of these objects are visible even in 5435-band, indicating that 
they reside at z < 4. 

On the other hand, the IR parts of their SEDs seem to be consistent with those of 
evolved stellar populations. Indeed, a SED template of local E/SO (e.g., Coleman, Wu & 
Weedman 1980) redshifted to 2; ~ 1-3 can reasonably fit most of the objects in the IR 
region. Such an empirical template, however, has a major drawback that its age (> 11 
Gyr) is much older than the age of the universe at the inferred redshifts. To overcome this 
problem, we explored the stellar population synthesis models of Bruzual & Chariot (2003), 
which provides the flexibility of adjusting the parameters of a model galaxy, especially its age 
and star-formation history. We used a Chabrier IMF (Chabrier 2003) and solar metallicity. 
As an extremely red color can also be produced by dust, wc also considered the effect of dust 
using the extinction laws of Calzetti et al. (2000; at A < 2000A) and Fitzpatrick (1999; at 
A > 2OOOA). 

For most of the lEROs, we find that the shape of their SEDs in the 1 to 8 fim wavelength 
range is best explained by the presence of a well-evolved stellar component with ages in the 
range of 1.5-3.5 Gyr, observed at redshifts 1.6 < z < 2.9. The gentle curvature of the 
infrared SED is like that expected from older stars, whose fi^ flux density peaks around 
1.6/im in the rest frame (e.g., Simpson & Eisenhardt 1999). It is not well matched by 
models dominated by heavily dust-obscured, young starlight (Fig. 2, top panel). Although 
1.5-3.5 Gyr-old stars formed in an instantaneous burst (or Simple Stellar Population; SSP) 
match the infrared SED well, such a model significantly under-predicts the observed optical 
(rest-frame ultraviolet) fluxes of these objects. We find that this can be solved by adding 
a weaker component of younger stars, which we model as a secondary instantaneous burst 
with an age of ~ 0.1 Gyr (Fig. 2, bottom panel). This two-component scheme can explain 
the majority of the lEROs if the composite spectra are put at the right redshifts. The 
best fit for each object requires fine-tuning the ages of both components, but this level of 
adjustment is not well motivated here, since other variables may also be important, including 
more complex past star formation histories, metallicity, and dust extinction. For the sake of 
simplicity, we choose to fit all 17 lEROs by the combination of a 2.5 and a 0.1 Gyr-old bursts, 
leaving only the redshift and the burst amplitudes as adjustable parameters. We emphasize 
that this two-component model is a highly simplified approximation to the actual physical 
processes and should not be over-interpreted. In particular, this model does not exclude the 
possibility that dust reddening may play a role in the galaxy colors, but simply emphasizes 
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that the SED shape strongly favors the presence of a dominant older stellar population, 
formed at a significantly higher redshift than that at which the galaxy is observed. The 
young, secondary component represents weaker star formation that was on-going or recently 
completed at the observed redshift. The two-component model generally provides a much 
better fit than does the young, dusty model. As an example, the goodncss-of-fit for any of the 
top panel models in Fig. 2 is 3-4 x worse than the two-component fit shown in the bottom. 
The GOODS MIPS 24/im observations will further help distinguish if any of these lEROs are 
dusty. For instance, if they have only an insignificant amount of dust, our two-component 
model predicts that the brightest lEROs have 24//m flux densities only at ~ 2/xJy level and 
thus will not be seen even by the deep GOODS 24/xm observations. On the other hand, if 
the lEROs are very dusty, their 24/xm flux will be much higher and thus could possibly be 
detectable. 

We also note that there is a degeneracy between the photometric redshift obtained 
during the flt and the ages of the components (especially the primary components) . Choosing 
an instantaneous burst older than 2.5 Gyr as the primary component does not change Zp 
significantly (at a level of A2; ~ 0.1), but using a burst younger than 2.5 Gyr does tend to 
increase the derived Zp^ with /\z ~ 0.4 in extreme cases. However, the overall Zp values are 
not likely much different from the current fits, since the plausible templates can neither be 
much older than 2.5 Gyr because of the constraint from the age of the universe, nor much 
younger than 1.5 Gyr because of the shape of the SEDs. 

Fig. 3 shows the SEDs of all 17 lEROs (filled boxes), along with the fitted two- 
component models (open boxes connected by dashed lines). Four of these objects have 
very uncertain photometric redshifts. Object #10 (the faintest in the sample) is only signif- 
icantly detected in four of the bands considered here, and thus can be fitted by templates 
at a variety of redshifts. Object #1 and #2 do not have optical detections and can be 
fitted by templates over a large redshift range. They are noteworthy, however, as their esti- 
mated redshifts (3.6 and 3.4, respectively, albeit very uncertain), are the highest among the 
lEROs. The SED for object ^^5 is different from that of any other object in the sample, with 
steeply increasing flux density from 3.6 to 8;um, and yields the worst fit from our chosen 
models. Therefore, these four objects are excluded from most of the discussion below. The 
remaining 13 objects have Zp values ranging from 1.6 to 2.9, with the median at Zp — 2.4. 
The adopted model template is self-consistent with the redshift lower bound ol z — 1.6, 
because it gives r(3.6/2;) < 16 at 2; < 1.6, which does not meet our color selection criterion 
even considering systematic errors in photometry. We refer to the sample of these 13 ob- 
jects as the reflned sample. Their Zp values are listed in Table 2. This table also lists the 
corresponding rest-frame Kg absolute magnitudes of these objects, which are derived using 
the 3.6/xm measurements and the K-corrections based on the fltted model SED templates. 
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The masses of the primary, older stellar components in the refined sample range from 0.1 
to 1.6x10^^Mq (with a median value of 5x10^°Mq), which are 0.1 to 2x the characteristic 
M* = 8.3 X IO^^Mq at z = from Cole et al. (2001; rescaled to account for the choice of the 
Chabrier IMF that we adopt here). Using a younger template, for example, a 1.5 Gyr-old 
template, will decrease the mass estimate by no more than a factor of two. The masses 
of the secondary bursts are 2-3 dex less. The median estimated masses of these objects is 
~ 9x larger than the typical "best fit" stellar masses estimated by Papovich et al. (2001) for 
L* Lyman break galaxies (LBGs) at z ~ 3 (~ 6 x IO^Mq, again adjusting for the choice of 
IMF) and ~ 3x larger than the LBG masses derived by those authors using their "maximum 
M/L" models. 

Two of these lEROs, #7 and #9, are detected in X-ray by Chandra (Giacconi et al. 
2001; XID 515 and 605, respectively). The X-ray hardness ratios of both sources are very 
high, indicative of heavily absorbed emission. Using the derived Zp and assuming an intrinsic 
X-ray power-law slope of slope of F = 1.8, we estimate the obscuring column densities 
and the unobscured X-ray luminosities (Lx) in rest-frame 0.5-7kcV. Source ^7 has Lx 
~ 1.1 X 10'^'^ erg s~^ with log(iV//) ~ 23.5; for 7^9, which is detected solely in the hard 
2-7 keV band, we obtain Lx > 2 x 10"''^ erg s~^, assuming a la lower limit on intrinsic 
absorption log(A'^^) > 23.3. With X-ray luminosities far exceeding that of starburst galaxies 
{Lx ^ 10''^ erg s"^), these sources are hkely obscured (Type 2) AGN, but are below the QSO 
regime {Lx ^ 10^^'^ erg s"^; see Norman et al. 2002). We estimate their observer- frame 
3.6/xm AGN emission by extrapolating from the unobscured Lx assuming a fiat uFi, (e.g., 
Elvis et al. 1994). Even assuming neghgible rest- frame NIR extinction, the predicted 3.6//m 
contribution from the AGN would be 0.575 and 0.09 fiJy for ^7 and ^9, respectively, or 
11% and 1.5% of the measured lERO flux densities. We therefore conclude that the AGN 
contribution to their optical-infrared SED is safely ignored and the fitting described above 
is not affected by the AGN. 



4. Discussion 

It is instructive to consider the relation between lEROs and other classes of faint, red 
galaxies, such as conventional EROs. From the Kg magnitudes and the interpolated R 
magnitudes listed in Table 1, one can see that nearly all lEROs do satisfy the conventional 
ERO criterion of {R - K)ab > 3.35 mag (or {R — K)vega > 5 mag). Arc lEROs just EROs, 
but selected at redder wavelengths? The answer does not seem to be that simple, because we 
find that most bright EROs do not meet the lERO color criterion. Moustakas et al. (2004) 
present a sample of EROs in the CDF-S with Kg < 22 mag, 16 of which are in the HUDF 
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and match IRAC sources. However, none of these 16 EROs is selected as an lERO, either 

because they do not pass the color criterion or because their IRAC photometry is uncertain 
due to blending. As all but one of our lEROs have Kg > 22 mag, it is important to consider 
fainter ERO samples. We use interpolated i?-band magnitudes and the ISAAC Kg-hand 
data to select a sample of 27 fainter EROs at 22 < Kg < 24 mag. Two-thirds (18/27) of 
them match the initial list of 75 lERO candidates, including sources that we the rejected due 
to blending issues. Among these 18 fainter EROs, 13 match our ZERO sample, i.e., roughly 
50% (13/27) of this fainter sample of EROs quahfy as lEROs. However, another 33% (9/27) 
do not match any ZERO candidates. 

Therefore, while most bright {Kg < 22 mag) EROs are not lEROs, there is a substantial 
(but not complete) overlap between the two selection criteria at fainter magnitudes. A 
consistent interpretation would be that the ZERO color criterion picks up the fainter, higher 
redshift equivalents of conventional EROs. The photometric redshifts of our lEROs are at 
1.6 < z < 2.9 (possibly extending to ^ < 3.6), with the median at 2; ~ 2.4. By contrast, 
redshifts (photometric and spectroscopic) for the Moustakas ct al. ERO sample peak at 2; ~ 1 
(only a minority at 1.6 < z < 2.5), with the median at z ~ 1.2. Furthermore, the (^850, 
3.6/xm) filter pair samples roughly the same rest-frame wavelength range at 2; ~ 2.4 as does 
the {R, Kg) filter pair at ;2 ~ 1.2, and T{3.6/z) > 20 corresponds to (2^850 "'"^3.6 ) > 3.25 mag, 
i.e., the ZERO criterion selects roughly the the same rest-frame feature at higher redshifts 
as the ERO criterion does at lower redshifts. 

Several authors (e.g., Totani et al. 2001; Franx et al. 2003; Saracco et al. 2004; Chen & 
Marzke 2004) have used near-infrared colors to identify high redshift galaxy candidates. All of 
our ZEROs with sufficiently deep J- and if -band photometry satisfy the (J — K)vega > 2.3 
mag criterion used by FYanx et al., and have a similar surface density on the sky. Five 
ZEROs also satisfy the (J — K)vega > 3 mag criterion used by Totani et al. and Saracco et 
al. Most of these objects populate the more distant end of our photometric redshift range 
(2.7 < Zp < 3.6). Chen & Marzke (2004) report nine faint, red galaxies selected from the 
HUDF ACS and NICMOS data, with (2775 -i/ieo) > 2 mag and Zp > 2.5. Four of our lEROs 
(7^1, 2, 8 and 9) are also in the Chen & Marzke list. Most of these authors have suggested 
that these red, IR-selected galaxies have relatively massive, old stellar populations at z > 2, 
although Totani et al. interpret them as highly obscured, star-bursting proto-ellipticals. By 
extending photometry to 8/xm with ZRAC, we also find that the SED shapes for most of the 
ZEROs favor the presence of dominant, old stellar populations. 

Zf most of these lEROs are indeed at z 1.6-2.9, their IR luminosities indicate that 
they are at the bright end of the luminosity function. The rest-frame i^T^-band absolute 
magnitudes of the objects in the refined sample range from Mab{Ks) — —21.3 to —24.2 
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mag, while the median is at —22.9 mag. Most of these objects are significantly brighter 
than the near-IR M* value in the local universe. For example, the 2MASS /r^-band LF of 
Kochanek et al. (2001) gives M*{Ks) = —23.53 + 51og(/i) for early type galaxies, which 
corresponds to M\q{Ks) — —22 A3 mag. 

A detailed treatment of the space density of lEROs is beyond the scope of this paper, 
and requires a more extensive analysis of source detection efficiency and incompleteness in 
the crowded IRAC images, which is presently underway. Here, we limit ourselves to some 
simple comparisons to the number density of local, early-type galaxies, which we regard 
as likely descendants of the lEROs. We assume the volume over which lEROs can be 
observed extends from z = 1.6, consistent with the lowest redshift at which an unreddened 
old stellar population should meet the lERO color criterion, to z — 2.9, the highest redshift 
for an object in the refined sample. The least luminous lERO (#12) in the refined sample 
has Mab{Ks) — —21.39 mag at Zp — 1.9. If we assume a roughly flat SED (in f^, vs. 
A) in the IRAC bands (as is typically observed), and redshift it to ^ = 2.9, this galaxy 
would have m(3.6;um) = 24.1 mag, and would be easily detected with a formal point source 
S/N~38. The actual detection likelihood will be smaller due to the blending, which affects 
our catalogs over a wide range of flux. We neglect this here, and naively assume that the 
refined sample represents a lower bound to a complete and volume-limited set of objects 
with MABiKs) < -21.3 mag at 1.6 < ^ < 2.9. 

Within this volume, the refined lERO sample with Mas < —21.3 mag has a space 
density of 3.2 x 10^^ Mpc~^, of which the six most luminous objects (with Mab{Ks) < 
—23.0 mag; hereafter "the luminous subsamplc") contribute 1.3 x 10~"^ Mpc^'^. The actual 
number density of such objects is certainly larger, due to catalog incompleteness and the 
conservative screening process during the selection. We compare these number densities to 
those calculated by integrating the local X^-band LF for early-type galaxies (Kochanek et 
al. 2001). For the (unphysical) case of no evolution (NE), the local LF predicts 12.4 and 
1.3 X 10-^ objects per Mpc^ for Mab{Ks) < —21.3 and < —23.0 mag, respectively. The 
(likely incomplete) lERO sample therefore makes up 25% of the number density of early-type 
galaxies with comparable luminosities in the local universe, while the luminous subsamplc 
has virtually the same space density as for comparably luminous local ellipticals. We might 
assume instead that the lERO population undergoes pure luminosity evolution (PLE) down 
to the present day. A simple, single-age stellar population fades by 1.0 mag in the rest-frame 
X-band from an age of 2.5 to 13 Gyr (Bruzual & Chariot 2003). In this scenario, the entire 
refined sample and the luminous subsample make up 14% and 20% of the corresponding 
local number densities of early-type galaxies with Mab{.Ks) < —20.3 and < —22.0 mag, 
respectively. We note that the galaxies in the luminous subsample are found entirely at the 
upper end of the redshift range, 2.4 < z < 2.9, and if this volume were used instead, their 
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space densities would be 2.55 x larger, making up 51% of the local population under the 
PLE model to Mab{Ks) < -22.0 mag. 

The HUDF samples a small cosmic volume, and we may expect the most luminous 
and massive galaxies to cluster strongly, potentially leading to large field-to-field variance 
in the densities. We must therefore regard these simple estimates with caution, pending 
analysis of the full GOODS area, with proper treatment of sample incompleteness (although 
we note that only the HUDF has sufficiently deep optical imaging data to enable firm color 
constraints on lEROs down to the flux limit of the Spitzer IRAC images). Nevertheless, 
we may say that the lEROs appear to be a signiflcant and numerous population of objects, 
most naturally interpreted as galaxies with relatively old and massive stellar populations at 
z > 1.6. 

To summarize, using the first epoch of the GOODS SpitzerLegacj Program observations 
in the CDF-S, we have identified 17 well-isolated objects in the HUDF that are significantly 
detected by IRAC but are very faint in the ACS images. Their SEDs are best explained by 
the presence of an old (~ 1.5 to 2.5 Gyr) stellar population in galaxies at 1.6 < 2; < 2.9. A few 
of the objects may have higher redshifts, but this is yet uncertain. The old stars dominate the 
infrared light, with a median rest-frame X^-band luminosity that is 1.5 x that of present-day 
L* early-type galaxies, and stellar masses ~ 0.1 to 1.6 x IO^^Mq for a Chabrier IMF. They 
are substantially more massive than typical Lyman break galaxies at similar redshifts. A 
much smaller component of recent star formation is needed to explain the optical (rest-frame 
UV) portion of the SED. The lEROs arc likely the highcr-redshift analogs of conventional 
EROs, selected via similar color criteria applied at longer wavelengths. The GOODS IRAC 
observations contribute to the mounting evidence for a significant population of red, evolved 
galaxies at high redshifts. In a simple PLE scenario, the lEROs may be direct progenitors 
for at least 14 to 51% of the local population of massive, early-type galaxies. 

Support for this work, part of the Spitzer Space Telescope Legacy Science Program, 
was provided by NASA through Contract Number 1224666 issued by the Jet Propulsion 
Laboratory, California Institute of Technology under NASA contract 1407. 
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Fig. 1. — As an example of IRAC-selected extremely red objects, 11-band image cut-outs 
of object #4 in Table 1 are shown here. Images are 12 xl2 with N up and E to the left. 
The optical (top panel) BViz images are from the HUDF ACS campaign. The near-infrared 
(middle panel) JH images are from the HUDF NICMOS Treasury Program, while the Ks 
image is from the data obtained at the VLT/ISAAC as part of the GOODS ground-based 
supporting observations. The IRAC images (bottom panel) are described in this paper. The 
location of this source, as derived based on the IRAC 3.6;um image, is illustrated with a l" 
radius circle in the top and the middle panels. 
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Fig. 2. — {Top) The SEDs of the lEROs cannot be easily explained without a well-evolved 
stellar population. The case for #8 is demonstrated here. The observed SED is shown 
as filled boxes while the template is shown as open boxes connected by dashed lines. The 
template is a 0.1 Gyr-old star-bursting galaxy, reddened by dust extinction of Ay = 1.0, 2.0 
and 5.0 mag for the left, middle and right panels, respectively. The corresponding best-fit 
Zp values are also labeled. (Bottom) The SEDs of most of the lEROs can be satisfactorily 
explained by a two-component model, and the fitting is shown for the same object ^8. 
The primary and the secondary components of the model are shown in the left and the 
middle panels, respectively. The primary component is a massive, old stellar population 
approximated by a 2.5 Gyr-old instantaneous burst, while the secondary component is a 
weak, on-going star-forming population approximated by a 0.1 Gyr-old instantaneous burst. 
The right panel shows the composition of these two components. The best-fit photometric 
redshift is 2;„ = 2.9. 
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Fig. 3. — For simplicity, a composite model of a 2.5 Gyr-old burst and a 0.1 Gyr-old burst 
is used to fit all the 17 objects. Legends are the same as in Fig. 2. Flux density upper 
limits, when available, are shown as crosses. For clarity, the SEDs of these objects have been 
brought to the same scale. The objects in the refined sample (see text) have Zp range from 
1.6 to 2.9, with the median at 2.4. The median rest-frame i^^-band absolute magnitude thus 
derived is -22.9 mag, which is significantly brighter than the M\^{Ks) value of the local 
elliptical galaxy LF. 



Table 1. Photometric properties of lEROs selected in the HUDF 



ID 




RA&DBC(J2000) 




S435 


Veoe 


(R) 




»775 






2850 




'.^110 


^160 


Ks 


m(3. 6^m) 


m(4.5^m) 


m(5.8^(m) 


m(8.0^tm) 


r(3.6/z) 


1 


3 


32 


42 


86 


-27 


48 


09 


29 


> 29.8 


> 30.2 






> 30.1 


> 29.5 


27. Olio. 10 


24.63±0.02 


23.63±0.09 


21 


99±0 


01 


21 


86±0 


02 


21.35±0.07 


21.39±0.09 


> 1009 


2 


3 


32 


38 


74 


-27 


48 


39 


96 


> 29.8 


> 30.2 






> 30.1 


> 29.4 


26.58±0.05 


24.68±0.01 


24.28±0.24 


22 


10±0 


01 


21 


88±0 


01 


21.67±0.05 


21.78±0.07 


> 831 


3 


3 


32 


43 


62 


-27 


46 


38 


96 


30.14±0.48 


28.14±0.06 


(28 


16) 


28 


18±0 


07 


27 


91±0 


09 






23.66±0.18 


22 


28±0 


01 


21 


79±0 


01 


21.38±0.05 


21.41±0.05 


178.8 


4 


3 


32 


41 


74 


-27 


48 


24 


84 


> 31.1 


30.22±0.36 


(29 


54) 


28 


88±0 


12 


27 


96±0 


09 


25.88±0.03 


24.60±0.01 


23.82±0.11 


22 


40±0 


01 


22 


26±0 


02 


22.30±0.12 


22.69±0.20 


167.3 


5 


3 


32 


43 


66 


-27 


48 


50 


69 


29.71±0.23 


29.33±0.11 


(29 


32) 


29 


29±0 


12 


28 


76±0 


13 






> 25.0 


23 


28±0 


03 


22 


92±0 


04 


22.37±0.11 


21.19±0.05 


156.6 


6 


3 


32 


33 


25 


-27 


47 


52 


19 


29.75±0.50 


28.96±0.17 


(28 


60) 


28 


15±0 


09 


27 


20±0 


06 


26.14±0.05 


24.88±0.02 


24.17±0.18 


22 


5()±0 


02 


22 


69±0 


04 


22.77±0.24 


23.22±0.40 


75.4 


7 


3 


32 


32 


16 


-27 


46 


51 


60 


28.72±0.16 


27.36±0.03 


(27 


19) 


26 


94 ±0 


03 


26 


71±0 


04 






23.73±0.15 


22 


15±0 


01 


21 


93±() 


02 


21.63±0.07 


21.61±0.07 


67.1 


8 


3 


32 


35 


06 


-27 


46 


47 


46 


27.94±0.11 


26.67±0.02 


(26 


51) 


26 


26±0 


02 


25 


86±0 


02 


25.11±0.04 


23.42±0.01 


22.82±0.06 


21 


44±0 


01 


21 


37±0 


01 


21.01±0.05 


21.14±0.06 


58.6 


9 


3 


32 


39 


16 


-27 


48 


32 


44 


> 30.8 


27.52±0.04 


(27 


01) 


26 


44±0 


02 


26 


12±0 


02 


25.45±0.03 


23.36±0.01 


22.87±0.06 


21 


93±0 


01 


21 


96±0 


02 


21.74±0.10 


22.10±0.16 


47.6 


10 


3 


32 


37 


16 


-27 


48 


23 


54 


> 31.4 


> 30.2 


(30 


38) 


29 


26±0 


12 


28 


49±0 


10 


> 26.5 


> 26.5 


> 26.5 


24 


38±0 


09 


24 


62±0 


21 


> 23.9 


> 23.8 


44.3 


11 


3 


32 


30 


26 


-27 


47 


58 


24 


29.81±0.50 


28.50±0.10 


(28 


06) 


27 


55±0 


05 


26 


96±0 


05 






23.87±0.14 


22 


96±0 


01 


22 


83±0 


02 


22.86±0.10 


22.91±0.12 


39.8 


12 


3 


32 


29 


82 


-27 


47 


43 


30 


30.11±0.45 


27.87±0.04 


(27 


53) 


27 


10±0 


02 


26 


98±0 


03 






24.18±0.14 


23 


30±0 


03 


23 


31±0 


05 


23.24±0.27 


23.71±0.44 


29.4 


13 


3 


32 


3!) 


11 


-27 


47 


51 


61 


28.87±0.16 


27.11±0.02 


(26 


79) 


26 


37±0 


01 


25 


57±0 


01 


24.40±0.01 


23.23±().()1 


22.64±0.05 


21 


98±0 


01 


22 


()3±0 


02 


21.87±0.08 


22.68±0.18 


27.2 


14 


3 


32 


48 




-27 


47 


07 


58 


2!).32±0.45 


27.12±().()4 


(26 


35) 


25 


63±0 


01 


24 


73 ±0 


01 






21.89±0.04 


21 


30±() 


01 


21 


27±() 


01 


21.7!)±0.10 


22.47±0.20 


23.6 


15 


3 


32 


38 


76 


-27 


48 


27 


07 


28.54±0.18 


27.21±0.04 


(27 


12) 


26 


97±0 


03 


26 


84±0 


05 


26.27±0.04 


24.79±0.02 


24.56±0.23 


23 


41±0 


04 


23 


34±0 


06 


23.32±0.31 


23.12±0.29 


23.6 


16 


3 


32 


36 


71 


-27 


46 


38 


96 


27.26±0.04 


26.68±0.02 


(26 


40) 


26 


02±0 


01 


25 


54±0 


01 


24.94±0.03 


24.10±0.02 


23.34±0.08 


22 


14±0 


01 


22 


10±0 


01 


22.00±0.08 


22.08±0.08 


23.0 


17 


3 


32 


33 


67 


-27 


47 


51 


04 


28.45±0.19 


26.76±0.03 


(26 


27) 


25 


72±0 


01 


24 


97±0 


01 


23.56±0.01 


22.67±0.00 


22.39±0.04 


21 


71±0 


01 


21 


76±0 


01 


21.95±0.07 


22.24±0.10 


20.1 



"'^'The magnitudes are in AB system, which are related to flux density fu (in erg s ^ cm" -2 Hz-^) by m. = -2.5 X Ig(fu) - 48.60. 

^'The reported photometric errors of the IRAC bands reflect the random errors only. Typical systematic errors in these bands are at < 0.1 mag level. ' 

The JiiQ-band and H^QQ-hand magnitudes are adapted from the catalog released with the NICMOS HUDF data, which are also "MAG_AUTO" magnitudes extracted by SExtractor. Objects 5, 7 11, 12, and 13 

are outside of the NICMOS HUDF fleld; # 10 is within the NICMOS fleld but is not detected in either band. | 

"^'The JCg-band magnitudes are based on the deep ISAAC images obtained at the VLT, which are "MAG_AUTO" magnitudes extracted by SExtractor using the updated zeropoints. 

^'The (-R) values are not measured but obtained by interpolating from Vgo6 ^nd 1775 to 6500A, and are listed here only for comparing against the conventional (R-K) ERO deflnition. 
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Table 2. Photometric redshifts and absolute magnitudes of lEROs 



ID 


m(3.6/im) 




M 


K-correction 


M{K,) 


1* 


21.99 


3.6 








2* 


22.10 


3.4 








3 


22.28 


2.9 


-23.22 


-0.17 


-23.39 


4 


22.40 


2.7 


-22.96 


-0.10 


-23.06 


5* 


23.28 


2.8 








6 


22.50 


2.3 


-22.56 


0.00 


-22.56 


7 


22.15 


2.7 


-23.22 


-0.10 


-23.32 


8 


21.44 


2.9 


-24.06 


-0.17 


-24.23 


9 


21.93 


2.8 


-23.51 


-0.14 


-23.65 


10* 


24.38 


2.1 








11 


22.96 


2.4 


-22.19 


-0.02 


-22.21 


12 


23.30 


1.9 


-21.39 


0.06 


-21.33 


13 


21.98 


1.9 


-22.71 


0.06 


-22.65 


14 


21.30 


1.6 


-23.05 


0.11 


-22.94 


15 


23.41 


2.7 


-21.96 


-0.10 


-22.06 


16 


22.14 


2.4 


-23.01 


-0.02 


-23.03 


17 


21.71 


1.6 


-22.64 


0.11 


-22.53 



The four objects marked with asterisks have very uncertain 
photometric redshifts and are not included in the refined sample. 

•^■Thc M values arc the absolute magnitudes at the rest-frame 
wavelength that correspond to the observer-frame 3.6/im, while 
M{Ks) values are the absolute magnitudes at the rest-frame Ks- 
band obtained by adding the K-corrections to the M. The K- 
corrections are calculated using the fitting templates as shown in 
Fig. 3. 



